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To advance quantum information science a constant 
pursuit is the search for physical systems that meet the 
stringent requirements for creating and preserving quan¬ 
tum entanglement. In atomic physics, robust two-qubit 
entanglement is typically achieved by strong, long-range 
interactions in the form of Coulomb interactions between 
ions or dipolar interactions between Rydberg atoms [1-4] . 
While these interactions allow fast gates, atoms subject 
to these interactions must overcome the associated cou¬ 
pling to the environment and cross-talk among qubits [5- 
8]. Local interactions, such as those requiring significant 
wavefunction overlap, can alleviate these detrimental ef¬ 
fects yet present a new challenge: To distribute entan¬ 
glement, qubits must be transported, merged for interac¬ 
tion, and then isolated for storage and subsequent opera¬ 
tions. Here we show how, via a mobile optical tweezer, it 
is possible to prepare and locally entangle two ultracold 
neutral atoms, and then separate them while preserving 
their entanglement [ - 1]. While ground-state neutral 
atom experiments have measured dynamics consistent 
with spin entanglement [10, 12, 13], and detected en¬ 
tanglement with macroscopic observables [14, 15], we are 
now able to demonstrate position-resolved two-particle 
coherence via application of a local gradient and par¬ 
ity measurements [1]; this new entanglement-verification 
protocol could be applied to arbitrary spin-entangled 
states of spatially-separated atoms [16, 17]. The lo¬ 
cal entangling operation is achieved via ultracold spin- 
exchange interactions [9-11], and quantum tunneling is 
used to combine and separate atoms. Our toolset pro¬ 
vides a framework for dynamically entangling remote 
qubits via local operations within a large-scale quantum 
register. 

Internal spin states of particles provide a robust and 
long-lived storage for quantum information. While engi¬ 
neered spin-dependent interactions can realize entangling 
gates between spins, they also predispose the system to 
strong environmental coupling and decoherence [9] . Spin- 
exchange interactions, which arise from a combination of 
quantum statistics and spin-independent forces, afford 
a promising alternate route to entanglement generation, 
and have been explored with both electrons in quantum 
dots and atoms [9-11, 18-20]. When two particles inter¬ 
act, their interaction energy depends on the spatial sym¬ 
metry of the two-particle wave function. If the particles 
have spin but are otherwise identical, the symmetry of 
the two-particle spin state directly determines the spatial 
symmetry of their wave function: For repulsive interac¬ 
tions, two bosons (fermions) of opposite spin in a triplet 


configuration experience enhanced (suppressed) interac¬ 
tions, while the converse occurs for the singlet spin state 
(Figure lb). By preparing a superposition of the triplet 
and singlet, dynamical quantum beats result in the ex¬ 
change of spin between the particles. In our experiment 
we prepare atoms of opposing spin in the lowest two mo¬ 
tional states (e and g ) of an optical tweezer potential that 
we represent as | which results in equal population 

of the singlet and triplet spatial wave functions ^(ri, r 2 ) 
and Vr(iT, r 2 ). The difference in the contact interaction 
energy between these states yields spin-exchange dynam¬ 
ics at a rate J ex , which depends on the s-wave scattering 
length and two-particle density (see supplementary ma¬ 
terials). Specifically, it gives rise to the dynamics 

IV’(i)) = \ te,ig) cos(J ex t/2h)+i\ie,tg) sin(J ex t/2h). (1) 

This evolution is associated with the effective spin- 
dependent Hamiltonian H = J ex S e -S g , where S e and S g 
are the spin-operators for the respective motional states. 
A spin-entangled state (| '\ e ^g) + i\ leitg))/ can be 
created by allowing the state to evolve for an exchange 
time of 7iA/2J ex . 

We schematically represent the experiment in Fig¬ 
ure la, in which separated optical tweezers on the left 
(L) and right (R) each containing a single atom are 
dynamically reconfigured to produce spin-exchange dy¬ 
namics and non-local entanglement. Atoms are com¬ 
bined in the right optical tweezer where spin exchange 
creates the entangled state (| t e,lg)R + A \ref\g) r)/ \/2. 
Importantly, in our experiments we convert this entan¬ 
glement into spatial-spin correlations by separating the 
atoms into two tweezers (Figure la) to yield a state 
(I T*>l| 1 9 )r + i\ U)l\ t^)i?)/V2. While verification of 
such entanglement is a standard tool in ion and Ryd¬ 
berg experiments [1, 3, 4, 21], spatially-resolved detec¬ 
tion of the entanglement present in interacting systems 
of ground-state neutral atoms is challenging. Theoreti¬ 
cal proposals have studied ways of detecting spatial [22] 
and spin [16] entanglement, and very recently experi¬ 
mental progress has been made using a quantum gas 
microscope [17]. We devised a protocol that yields ro¬ 
tations of a two-qubit entangled state on the associated 
Bloch sphere via a combination of a magnetic-field gradi¬ 
ent and global spin rotations (Figure lc), allowing detec¬ 
tion of basis-independent correlations for arbitrary entan- 
gled states of the form + e^|| g )L|t 9 }fl)- 

Our protocol is applicable to qubit pairs in a large quan¬ 
tum register [11], to interacting spins in a Bose-Hubbard 
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FIG. 1. Experimental overview, (a) Two particles, trapped in separate optical tweezers, are initiated in their motional ground 
state and in opposing spin states (left). Dynamical reconfiguration of optical tweezer traps merges the particles into the same 
optical tweezer in specific motional states, and spin exchange entangles the particles. The atoms are separated into different 
optical tweezers, and the entanglement between the particles is experimentally detected (right), (b) Exchange interaction. Two 
atoms of opposite spin are prepared in the ground (g) and excited (e) motional states of the x direction. For bosonic atoms the 
singlet channel anti-symmetrizes in the motional degree of freedom, yielding 'ips(^1^2) = ('0^(ri)'0 e (i’2 ) - i>e(r 1 )ipg{r 2 ))/V 2 . 
From one-dimensional spatial wave functions we see that the anti-symmetrization prevents the two particles at positions x\ and 
X 2 (in units of the oscillator length 7 * 0 ) from occupying the same position (dashed white line), which is true for all choices of 
y, z. This leads to zero interaction energy for the singlet channel and J ex for the triplet channel with associated wave function 
'tpT (r 1 , r 2 ) = (^g( ri)'0 e (r2) + ( 1*1 )i/j g (r2))/^/2. (c) Studying entanglement of the separated particles. On the Bloch sphere, 
the exchange interactions ideally result in a state pointing along the zb y-axis, but in general can point in any direction in the 
equatorial plane with associated coherence angle 4>. 


chain [13, 16], and to strongly interacting fermions fea¬ 
turing anti-ferromagnetic correlations [23-21 ]. 

The experiment begins by loading two thermal 87 Rb 
atoms into two separate optical tweezer potentials [26]. 
Each atom is then separately laser cooled to the 3D 
ground state via Raman-sideband cooling, leaving a 3D 
ground state fraction of 90(7)% [ , 27, 28]. The atoms 
are initialized in opposite spin states with a fidelity in the 
range of 80 — 91% via single-atom addressing (see sup¬ 
plementary materials). We then prepare the pure two- 
particle state | t e,lg)R o n a single tweezer by reducing 
the separation between the two optical tweezers, and us¬ 
ing our capability to control the tunnel coupling between 
wavelength-scale optical tweezer traps [24, 28]. As illus¬ 
trated in Figure 2a, the spacing between the two tweezers 
is decreased until tunneling is appreciable; however, un¬ 
like our previous work [28], we dynamically shift to an 
asymmetric configuration such that the ground state of 
one optical tweezer is near resonant with the first (ra¬ 
dial) excited state of the other optical tweezer. We then 
perform adiabatic passage by slowly tuning the relative 
tweezer depths linearly in time, which coherently trans¬ 


fers the left atom into the right optical tweezer over the 
12 ms duration of the ramp (see supplementary materi¬ 
als). 

After a desired evolution time in the presence of ex¬ 
change dynamics, the adiabatic passage is applied in re¬ 
verse, yielding a motional state mapping of the excited- 
state atom back into the ground state of the left well. We 
can then read out the spin populations to verify the ex¬ 
change oscillations by ejecting atoms in the |t) state from 
the optical tweezers, and imaging the atom occupancy in 
each of the two tweezers. With this procedure, we can 
ascertain what the spin and motional degrees of freedom 
of each of the two atoms were when they occupied the 
same optical tweezer. 

Exchange oscillations in our experiment are shown in 
Figure 2b, and show the expected anti-correlated behav¬ 
ior. For these data, we remove experiments in which 
imperfections in our spin preparation lead to the spins 
remaining aligned; such events yield a static contribu¬ 
tion to the signal. We observe undamped oscillations out 
to times as long as 100 ms, despite the fact that our 
single-particle coherence time between |j") and |^) is less 
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FIG. 2. Direct observation of spin-exchange dynamics between two atoms, (a) Preparation of motional state configuration and 
detection in the double-well potential, (b) Using post-selection on our spin preparation, we plot the likelihood to measure the 
state \^ g )L\ig)R (green) and \i 9 )L\tg)R (purple) as a function of time between the end of the first adiabatic passage and the 
beginning of the second, (c) Measured spin-exchange oscillation frequency as a function of optical tweezer depth. The dashed 
black line is the predicted exchange frequency from a parameter-free model of the potential (see supplementary materials); the 
shaded region indicates the effect of systematic uncertainties on this prediction. All error bars are the standard error in the 
measurement. 


than 1 ms due to magnetic-field fluctuations. This is an 
expected feature of the entangled states created by the 
exchange interaction: homogeneous magnetic-field fluc¬ 
tuations induce a global phase on the two-particle super¬ 
position and, as such, leave quantum measurements un¬ 
affected. Hence, the state occupies a so-called “decoher¬ 
ence free subspace” [10, 29]. We can also control the fre¬ 
quency of the spin oscillations by modifying the depth of 
the optical tweezer in which the exchange occurs, which 
tunes the two-particle density. To study this dependence, 
we prepare | t e^g)R and linearly increase the depth of 
the tweezer in 5 ms, allow evolution of the exchange, and 
then ramp back in reverse and perform the second adia¬ 
batic passage. We can model the three-dimensional non- 
separable potential of our optical tweezer trap [30], and 
find agreement (Figure 2c) between the calculated and 
measured spin-exchange frequency J ex /(27rfi) (see sup¬ 
plementary materials). 

The measurements presented thus far have shown cor¬ 
relations in single-particle spin states. However, to en¬ 


sure that future operations can retain and propagate 
quantum information, one must verify that the phase 
coherence within the entangled state is preserved upon 
separating the particles. The entanglement verification 
protocol for separated atoms is summarized in Figure 3a. 
For explanatory purposes, we first focus on the case 
when the particles are separated after an exchange time 
of t ent = nivh/2J e x where n is an odd integer. The 
entangled state after the second adiabatic passage is 
IV’i) = t)zj \)r ± *| I)l| T)ij), omitting from 

now on the ground-state (g) motional subscripts to sim¬ 
plify notation. The \'ip±) states correspond to the gray 
and orange Bloch vectors, respectively, in Figure 3b. We 
then apply a magnetic-field gradient that imposes a dif¬ 
ference, Sh, in the | t) \ 1} single-atom-transition 
energy between the left and right optical tweezer. By 
applying the gradient for a time t g , a transformation 
|V>±) ->• ± ie lSt s\l) L \ t)i?) is achieved. As a 

function of t g , the state rotates between the singlet (pink 
in Figure 3b) and triplet (blue) with frequency 5. We 
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FIG. 3. Detection of non-local entanglement, (a) Procedure for creating and detecting entanglement, (b) Highlight of four 
different Bloch vector orientations in the experiment. The gray and orange orientations correspond to the outcome of the 
spin-exchange dynamics (after the second adiabatic passage), which yield rotations in the y-z plane of the Bloch sphere. The 
pink and blue orientations correspond to the points of peak parity, and are accessed by applying a magnetic field gradient that 
rotates the states about the z- axis, (c) After creating |-0+), we plot the measured parity as a function of gradient time t g . The 
gray bar is the bound on parity oscillation contrast delineating separable and entangled states, which accounts for imperfect 
spin-preparation. The colored dashed lines are at times when the associated states indicated in (b) are created, (d) Measured 
parity as the exchange time is varied and correspondingly the atoms are entangled and unentangled. In the lower plot, we set 
t g such that it rotates to | T) and then measure the parity. The upper plot is the same experiment without the parity 
detection, i.e. the protocol of Figure 2. The dashed lines indicate times when the corresponding states indicated in (b) are 
produced. The error bars in the data plots are the standard error, and the pink swaths show the 95% confidence bands. 


then apply a global 7 t/ 2 pulse in the (|t), \i)} sub-space. 
This pulse maps the singlet back to itself, while it maps 
the triplet to a Bell state ^(| Z)l\ t)fl + I 1)l\ 1)r)- 
Therefore, by measuring the probability that the spins 
are aligned or anti-aligned as a function of t g , we can 
observe singlet-triplet oscillations whose amplitude char¬ 
acterizes the two-particle coherence. We quantify this 
probability with the parity n (t g ) = ZZj Pj(~ l)* 7 ', where 
Pj is the likelihood to measure j atoms in the spin-down 
state [] -4]. The parity is equivalently the projection of 
the Bloch vector in Figure 3b onto the x-axis prior to the 
7r/2-pulse, and hence the gradient is essential because, 
though entangled, the states \'ip±) (gray, orange) exhibit 
zero parity after application of a 7r/2-pulse. 

We demonstrate the outcome of the verification proto¬ 
col on the state |^+) in Figure 3c. We plot H(t g ) after the 
microwave 7r/2 pulse, and observe oscillations in the par¬ 
ity signal as the gradient time t g is scanned. The contrast 
of these oscillations is consistent with what is expected 


given the exchange oscillation contrast in Figure 2, and 
non-vanishing parity oscillation would certify entangle¬ 
ment in the ideal case of perfect spin preparation. How¬ 
ever, we have imperfect spin preparation, and the erro¬ 
neous spin populations outside the {IDzJt)#? 
manifold could lead to parity oscillations even in the ab¬ 
sence of entanglement. We have derived a condition on 
the parity oscillation contrast that is necessary and suffi¬ 
cient to certify entanglement and is the simplest way to 
see there is entanglement in our system (for a full deriva¬ 
tion of this condition and its relation to other entangle¬ 
ment metrics see supplementary materials). We relate 
the measured parity contrast, C, to the measured proba¬ 
bilities (ptt^pFi) that the spins are erroneously prepared 
in the same spin-state: If C > Cbnd = 4(P^P^) 1 / 2 , 
then the state is entangled. By directly measuring the 
spin populations (see supplementary materials) and their 
associated uncertainty, we ascertain Cbnd = 0.133(25) as 
indicated by the dashed lines in Figure 3c. The observed 
parity oscillation contrast C = 0.49(4) exceeds Cbnd by 
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more than 7cr, certifying the presence of entanglement 
in the final state of the separated spins. We verify en¬ 
tanglement without correcting the measured parity for 
experimental imperfections, such as single atom loss due 
to background collisions. 

While in Figure 3c we varied the parity detection pa¬ 
rameters via t g , in Figure 3d we measure the dependence 
of the parity on the exchange time at fixed t g , thereby 
observing oscillations as the exchange interactions peri¬ 
odically entangle and unentangle the two atoms. We fix 
t g in the parity detection such that the entangled state 
|i/q_) (gray lines in Figure 3b, d) is rotated to a peak in 
II, corresponding to the creation of the triplet (blue lines 
in Figure 3b, c). Because this t g amounts to a tt/2 ro¬ 
tation about the z-axis of the Bloch sphere, it will also 
rotate |^-) to the singlet, which corresponds to maxi¬ 
mal negative parity. In the lower panel of Figure 3d, 
we show how the parity measured under these condi¬ 
tions oscillates at the exchange frequency J ex /(27 tK) . For 
comparison, in the upper panel, we show the measured 
exchange oscillations (purple, green) without the parity 
detection. At the linear points of the exchange oscilla¬ 
tions, one expects maximal entanglement corresponding 
to states \f>+) (gray) and | ^_) (orange) and thus the ex¬ 
tremal points of the parity. At the minima and maxima 
of the exchange oscillations, the atoms are unentangled 


and the parity vanishes. 

In conclusion, we have demonstrated the entanglement 
of remote qubits using spin-exchange interactions and a 
general protocol for detecting spin entanglement in a di¬ 
versity of systems. In a large register, tuning the en¬ 
tanglement phase could be achieved by a far-detuned 
focused probe that changes the local effective magnetic 
field experienced by a single qubit, and the qubits can 
be arranged to allow the passage of the optical tweez¬ 
ers without perturbing the qubits [11, 31]. While in this 
work we focus on quantum information applications, our 
abilities to control spin and motion of individual neu¬ 
tral atoms will allow the study of intriguing microscopic 
models in condensed-matter physics, such as the Kondo 
Lattice model [32]. 
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Supplementary Materials 


I. STATE PREPARATION AND 
EXPERIMENTAL PROTOCOL 

The experiment begins with atoms loaded stochasti¬ 
cally from a magneto-optical trap into a pair of optical 
tweezers separated by 1.57 /im, each with an optical waist 
of 710 nm and a 23(1) MHz depth [26]. For all data pre¬ 
sented in the figures, we post-select on experiments in 
which each optical tweezer is loaded with a single atom, 
based upon atom population (0 or 1) measurements in 
each tweezer at the start of the experiment. We perform 
optical pumping and three-dimensional Raman sideband 
cooling as described in Refs. [27, 28], leaving each atom 
in the motional ground state of its respective tweezer 
with 0.90(7) fidelity, as determined via motional side¬ 
band spectroscopy, and in the |F, mp) = |2,2) = | t) 
spin state. To perform single-atom addressing, we subse¬ 
quently apply circularly-polarized light to the left tweezer 
to induce a relative shift between the tweezers of the 
|t) 11,1) = |t) transition [28, 3 ]. For the data in 
Figure 2b and Figure 3d, we apply a resonant square mi¬ 
crowave 7r-pulse to the non-light-shifted atom, yielding a 
preparation fidelity in the desired two-particle spin con¬ 
figuration of 0.835(10). For the data in Figure 3c, we 
improve this fidelity by using temporally shaped Gaus¬ 
sian pulses in order to minimize Cbnd via the suppression 


of off-resonant transitions. 

After motional and spin preparation, the position of 
the left optical tweezer is swept adiabatically in 10 ms to 
the other tweezer, realizing a two-atom spacing (Gaus¬ 
sian centroid separation) of 715 nm (854 nm). The 
tweezer depths are then reduced to 91(4) kHz in an adia¬ 
batic ramp [28] ; in this configuration, there is a measured 
resonant tunneling of J eg /(27 tK) = 165(6) Hz between the 
left-tweezer ground state and right-tweezer excited state. 

We then apply the first adiabatic passage (AP) through 
the tunneling resonance. Starting with the left optical 
tweezer tuned 2.2 kHz below the tunneling resonance, the 
left tweezer depth is swept linearly in 12 ms to symmet¬ 
rically above the resonance. The end of this ramp marks 
the nominal t = 0 of the exchange dynamics, and the 
resulting tweezer depths are static during the exchange 
time. After the exchange time, the AP is performed ex¬ 
actly in reverse, transferring the motional-excited atom 
back into its origin tweezer. The full AP procedure in 
both directions has a measured fidelity of 0.81(4) for the 
data in Figure 2b; we observe systematic fluctuations in 
this fidelity such that its value is 0.69(2) for the data in 
Figure 3c. 

At this point, if we are directly measuring the exchange 
dynamics as in Figure 2, the tweezers are ramped up in 
depth and separated for spin detection and imaging. If we 
are performing entanglement detection, instead directly 
we apply the magnetic-field gradient. After switching on 
the gradient there is a 3 ms delay for field settling. Then 
we hold the gradient on for variable time t g , after which 
the gradient is switched off. After another 3 ms delay to 
let the fields settle, we apply the microwave tt/2 pulse for 
parity read-out, after which the tweezers are increased in 
depth and separated for spin detection and imaging. 

The observed spin-exchange contrast can be influenced 
by a number of factors, including the single-particle 
ground-state fraction, single-particle loss, and the AP fi¬ 
delities. Spin preparation also affects the contrast, but 
as discussed in the text the data presented in Figure 2b 
are post-selected on the desired anti-aligned spin con¬ 
figuration (for the parity measurements, post-selection 
is not possible since we study all final spin configura¬ 
tions). Based on our measurements of experimental sys- 
tematics - the ARP fidelity of 0.81(4), the single-particle 
ground-state fraction of 0.90(7), and the single-particle 
survival probability of 0.963(7) - and combined with the 
spin-preparation post-selection, we would expect a spin- 
exchange contrast of 0.60(10), which exceeds the mea¬ 
sured value of 0.29(2) in Figure 2b. In the absence of 
the spin-preparation post-selection, we observe a spin- 
exchange contrast that is reduced by the spin-preparation 
fidelity of 0.835(10), as expected. 

For the characterization of the spin-populations, in a 
distinct experiment we measure the spin-down probabil¬ 
ity of the atom in each tweezer by applying resonant 
push-out light directly prior to when the first AP of the 
experiment would occur. Accounting for the separately 
measured single-atom survival probability of 0.963(7), we 
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measure {ptt^pU} = {0.071(14), 0.016(5)}, leading to 
pit _l_ ptl = l — (ptt p pH) . These numbers are used 
to calculate Cbnd and F. The measured parity data is 
not corrected for loss of single atoms due to background 
collisions, which degrades the entanglement generation. 


II. MODELING OF INTERACTIONS AND 
TUNNELING DURING STATE PREPARATION 

For two atoms of mass m in different motional wave 
functions ^(r) = (r|e) and V^( r ) = (r|g), the ex¬ 
change interaction energy J ex = 2 U eg , where U eg = 

47r ^ Qs f \ 2 p e (r)\ 2 \ 2 p g (r)\ 2 d 3 r is the contact interaction en¬ 
ergy proportional to a s , the atomic s- wave scattering 
length. To compute this integral we use efficient numeri¬ 
cal techniques to determine the spectrum and eigenstates 
of the 3D optical tweezer potential [30] using parameters 
from an independent experimental characterization [28]. 
This calculation, combined with experimental uncertain¬ 
ties, gives rise to the blue theoretical swath in Figure 2. 



FIG. SI. Adiabatic energy-eigenstates as a function of the 
double-well bias A in units of the ground-excited tunneling 
J eg . At large positive bias, the triplet and singlet eigenstates 
corresponding to two particles in the same well are split by 
Jex • 


To theoretically model the adiabatic passage proce¬ 
dure (AP), we consider a situation in which two par¬ 
ticles with opposite spins are initially localized in the 
ground motional states of the left and right wells, de¬ 
noted L g and R g , respectively. During the AP, the 
bias between wells A is tuned so that L g is near res¬ 
onance with an excited level of the right well, denoted 
R e . In what follows, we will measure the bias A with 
respect to this resonance position, i.e. the resonance oc¬ 
curs at A = 0. Further, we consider that the bias 
range is such that tunnel-couplings to all other motional 
states (e.g., L g —)> R g tunneling) are negligible, and so 
we can restrict ourselves to the set of motional states 
spanned by L g , R g , and R e and the single L g —>> R e 
tunneling resonance at A = 0. Because of the mixing 
of spin components due to the spin exchange interac¬ 
tion, it is most convenient to use the singlet-triplet basis 
{| S-,L g ,R g ), | S;R e ,R g ), \T;L g ,R g ), \T;R e ,R g )}, where 


S and T denote singlet and triplet states and the latter 
two labels are the motional states of the two particles. 
In this basis, the Hamiltonian is 

/ A -J eg 0 \ 

fr [ -Jeg 0 0 0 1 

0 0 A -Jeg I ’ 

\0 0 -Jeg 2U eg J 


where J eg is the tunneling amplitude for the process 
L g —>> R e , and the initial state at large negative A, \L g t 
, R g ^), is an equal weight superposition of \S]L g ,R g ) 
and \T]L g ,R g ). The AP process is described by sepa¬ 
rate tunneling avoided crossings in the singlet and triplet 
channels, as shown in Figure SI. Provided that the ramp¬ 
ing procedure is adiabatic with respect to these avoided 
crossings, it will transfer the initial state into an equal 
weight superposition of | S;R e ,R g ) and \T;R e ,R g ) at 
large positive A. This pair of eigenstates, which cor¬ 
respond to the two particles occupying the same well, 
have an asymptotic energy splitting of 2 U eg as | A| J eg . 
Note that the position of the bias resonance for the triplet 
channel is shifted by 2 U eg with respect to the resonance 
in the singlet channel. The energy offset does not affect 
the degree of adiabaticity of the ramping procedure and 
just gives rise to a phase shift between the singlet and 
triplet components. 


III. ENTANGLEMENT VERIFICATION BASED 
ON PARITY OSCILLATIONS 

Here we derive a criterion for verifying entanglement 
generated by spin-exchange interactions in a two-atom 
system. Our strategy will be to assume an unentangled 
(separable) density matrix, and from this assumption es¬ 
tablish a constraint on experimentally measurable quan¬ 
tities: the parity oscillation contrast (Fig. 3c) and the 
populations of different spin states. Experimental viola¬ 
tion of this constraint thus verifies entanglement. 

In the experiment, spin-exchange occurs between two 
atoms occupying a single optical tweezer, and then those 
atoms are separated into two tweezers. Because the adia¬ 
batic passages in the experiment are imperfect, the atoms 
may sometimes end up in the same tweezer after the at¬ 
tempted separation. For clarity of presentation, we first 
consider the idealized case of perfect adiabatic-passage 
fidelity; hence p in what follows describes states imme¬ 
diately after the second adiabatic passage in which one 
atom occupies each tweezer. At the end of this section we 
carefully consider the effects of adiabatic passage failure, 
and show that they do not affect our claims of entangle¬ 
ment. 

Because the measured spin-coherence time in the ex¬ 
periment is much less than the time between when we 
prepare the initial spin state and when we complete the 
adiabatic passages, p cannot have any coherences be¬ 
tween states with different total S z = S Z L + S Z R . Working 
in a basis that diagonalizes both S Z L and S R , {| t) l\ t 










)r, |t)L|l)i2, | >1)l | T)i ?5 |I)l|I)r}, the most general den- F = (T + |p|T + ). As shown in Ref. [1], p is entangled if 
sity matrix satisfying this condition can be written F > 1/2. In terms of Eq. (S2), the fidelity can be written 


(PAlAr) 0 o 

o PAlAr) £ 

0 £* P(IlAr) 

\o o o 


0 

0 

0 

P(IlAr) 


(S2) 


The populations PAlAr) and PAlAr) are, respec¬ 
tively, the total probabilities of having both atoms in 
the | t) state or both atoms in the | state, and are 
non-zero due to imperfect initial spin preparation. Be¬ 
cause these probabilities are conserved by the adiabatic 
passages and the spin-exchange, their measured values 
before the adiabatic passage, referred to as P^ and 
pU i n the manuscript, can safely be used in Eq. (S2): 
PAlAr) = P n and PAlAr) = P U • The parity is 
measured after first applying a magnetic-field gradient 
for a variable time t g , and then applying a 7r/2 microwave 
pulse, which transforms p p(t g ). After some algebra, 
the parity of p(t g ) can be written n(£ g ) = 2Re(e - ^ tg £), 
which oscillates as a function of t g with a contrast of 
C = 4|e|. 

Our goal is to derive a constraint on C in terms of 
the measured quantities pTf, P^, under the assumption 
that p is separable. If p were a product state pl ® Pr, 
where 


Pl{r) - 


t,t td \ 
Pl{r) Pl(r ) | 

1,1 I 5 

Pl(r) Pl(r ) / 


(S3) 


then 

N = I pP\ x \pP\ < (Pl^Pl 1 ) 1 ' 2 x (pPpP) 1 ^ (S4) 
= (P(tL,tfl)P(4-L,4.fl)) 1/2 = (P n P u ) 1/2 . (S5) 

By definition, a separable state can be written as a clas¬ 
sical mixture of product states, p = JP A ip° L (8) p J R , in 
which case 


M < £ N ^ (PpPf 1 ) 1/2 (S6) 

3 3 

< (E A j Pp) 1/2 {^Z xj Pp) 1/2 = (P n P u ) 1/2 . (S7) 

3 3 


F=*^[P(tL,l R ) + P(lL,U)] + j (S9) 

= \ + j -\{P n + P u )- (S10) 

Therefore P > 1/2 is equivalent to C > 2(pTf + P^), 
which agrees with Eq. (S8) when P ^ = P++. Equa¬ 
tion (S8) can also be derived by applying the Peres- 
Horodecki criterion to p, and thus it is actually both 
sufficient and necessary for entanglement [33, 34]. Thus, 
in contrast with the fidelity-based entanglement witness 
C > 2(PIT -f p44) ? the right-hand-side of Eq. (S8) is as 
small as possible for any PP and pTf, resulting in the 
greatest possible confidence in entanglement for a partic¬ 
ular measured contrast. We also note that the extent to 
which Eq. (S8) is satisfied, \{C — 4(pttpU^ 1 2 ) > q, is 
a direct measurement of the concurrence in the density 
matrix p [16]. 


A. Imperfect adiabatic passages 

As mentioned above, the adiabatic passages in the ex¬ 
periment are not perfect. For our purposes, we define 
success as any situation in which the atoms end up in 
different tweezers, which occurs either if both of the in¬ 
dividual adiabatic passages succeed or if they both fail. 
Conversely, we define failure as any situation in which 
both atoms end up in the same tweezer, which happens 
if one of the adiabatic passages is successful while the 
other is not. While the states resulting from failure have 
not been considered in deriving Eq. (S8), they do not con¬ 
tribute to the parity oscillation contrast, since two atoms 
in the same tweezer are not sensitive to a magnetic-field 
gradient. Therefore, we intuitively expect that imperfect 
adiabatic passage can only lower the measured contrast, 
such that Eq. (S8) still implies entanglement. 

This intuition can be formalized by introducing a pro¬ 
jector, iC, onto the states with one atom in each tweezer. 
Defining / as the success probability, we can then form 
projections of the true experimental density matrix, de¬ 
noted by p exp , into the subspaces defined by success or 
failure: 


Here, the second inequality is the constraint on |e:| derived 
for a product state in Eq. (S4) applied to each state in the 
classical mixture, and the first and third are triangle in¬ 
equalities. Using C = 4|e|, we therefore have guaranteed 
entanglement whenever 

C > 4(P n P u ) 1/2 . (S8) 

We can also verify entanglement by using an entangle¬ 
ment witness associated with the fidelity of p in the maxi- 
mally entangled state |$+) = (|t)LM.)fi+*M.)L|t)fl)/v / 2, 


Kp exp k = /ftmcc (Sll) 

(i - k) Pex p (i - k) = (i - f) Piaih (S12) 

Note that, by the choice of pre-factors, both p succ and 
Pfaii are properly normalized density matrices. Impor¬ 
tantly, Psucc has precisely the form given in Eq. (S2); we 
can therefore directly apply the arguments above, deduc¬ 
ing that Psucc is entangled whenever |^| > vPTfiPU. The 
only difference from before is in the relationship between 
|ej and the measured parity oscillation contrast. The ap¬ 
plication of a magnetic-field gradient and the subsequent 
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7r/2 pulse does not mix the two subspaces partitioned 
by K , and therefore the measured parity can be written 
n(tg) = /n s ucc(^) T (1 /)n fai i. Importantly, Ilfaii is 
independent of t g because two atoms in the same tweezer 
are not sensitive to a magnetic-field gradient, and there¬ 
fore the measured parity oscillation contrast is given sim¬ 
ply by C = / 4\e\. Thus the effect of imperfect adiabatic 
passages is that, in terms of the measured contrast, the 
condition |e:| > yfpttpU- now reads C/f> 4 y/pttp-U 
[c.f. Eq. (S8)]. The criterion used in the manuscript, 
C > 4\/pttp44, is therefore a conservative way to verify 
entanglement in p s ucc , since C/f is strictly larger than 

C. 

The fidelity can now be written F = (T + |p exp |T + ) = 
/(^+|Psucc|^+) = fF succ , where F succ is the fidelity of 
the (projected and normalized) density matrix p succ , and 
can be equivalently written as F succ = | + |er| — \{P^ + 


P^). Similar to before, entanglement of p succ is now 
guaranteed by F succ >1/2, or equivalently 4 |e:| > 2 (P^+ 
P u ). Using 4|er| = C/f, we can extract P succ from the 
measured contrast [C = 0.49(4)] and measured success 
probability [/ = 0.69(2)], obtaining F succ = 0.634(17). 
Note that the actual fidelity F is reduced from this value 
by the success probability /, and is in fact below 1/2. 
However, this is not inconsistent with entanglement in 
p exp , as the fidelity-based witness for p e x P , written in 
terms of the actual fidelity, is F > f / 2. 
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